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1 Introduction and main results 

A determining factor in the stability of the magnetic state of small ferromag- 
netic particles is the structure of the spectrum of their low-lying excitations. 
Stability against thermal (and quantum) fluctuations is a major concern 
when one is interested in increasing the density of information stored on 
magnetic hard disks. Higher density of information requires smaller mag- 
netic particles to store the bits. The smaller these particles get, the less sta- 
ble their magnetic state tends to be. It is also well-known that ferromagnets 
spontaneously form domains with different orientations of the magnetization. 
These two facts motivate us to study the excitation spectrum of finite size 
ferromagnets with a domain wall or interface. From examples, it is known 
that the presence of an interface, in general, has an effect on the low-lying 
excitation spectrum H, ^. 

We consider the spin 1/2 XXZ Heisenberg model on the three-dimensional 
lattice Z^. For any finite volume A C Z^, the Hamiltonian is given by 

Ha = ~ Yl ^-'(s^^s^^ + s^s^) + (1-1) 

\x-y\^l 

where A > 1 is the anisotropy. It will be convenient to work with the usual 
parametrization A = [q + q~^)/2, < g < 1. Note that in the limit A — »• cx) 
(g — > 0), one recovers the Ising model. The case A = 1 (g = 1) is the XXX 
Heisenberg model. 

It is well-known that this model has two ferromagnetically ordered trans- 
lation invariant ground states. What is less well-known is that there are also 
ground states describing an interface between two domains with opposite 
magnetization. The 100 interfaces are similar to the Dobrushin interfaces 
found in the Ising model. They exist for sufficiently small temperatures, as 
was recently proved in 0. Unlike the Ising model, the XXZ model also pos- 
sesses ground states with a rigid 111 interface at zero temperature p|. Its 
stability at positive temperatures is still an open problem. 

In this paper we are interested in estimating the low-lying excitations 
above the ground state with a 111 interface. It is easy to show that the 
excitation spectrum above the translation invariant ground states has a 
non-vanishing gap. In 0] it was proved that, in the corresponding two- 
dimensional model, the excitations above the 11 interface are gapless. By 
an extension of the methods in []TU[, Matsui |TT]] showed that the excitation 
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spectrum has to be gapless in all dimensions > 2. Here, we are interested 
in the nature of the low-lying excitations for the three-dimensional model, 
and in particular their dependence on size. We prove the following bound 
for the energy of an excitation localized in a finite domain of linear size R. 

Main Result: Excitations localized in A/j have a gap 'jji bounded by 

7/1 < 100^^^-^-^, for R>70. (1.2) 

where 5{q,v) is an exponent between and 1/2 that depends on the filling 
factor V of the interface plane (see explanation below), as well as the param- 
eter q. 



The meaning of this bound is the following. We consider the model 
in a finite volume A, with a fixed magnetization and boundary conditions 
that induce an interface. By perturbing the ground state in a cylindrical 
subvolume Ar, with circular cross-section of radius R, we then construct an 
orthogonal state with the same magnetization. The bound ( |1.2|) is an upper 
bound for the difference in energy of this state with respect to the ground 
state in the limit A y Z^. For finite volumes A, the same bound holds as 
long as A is substantially larger than R. When R and the finite volume are 
comparable in size, a similar bound holds but with a larger constant factor 
and additional error terms (see Section 4). 

The dependence on q of the bound ( |1.2| ) has some interesting features, 
which we explain next. First, in the limit g — 1, the bound diverges. This 
means that our Ansatz for the excitations of the 111 interface does not work 
for the isotropic model. This is not surprising as the isotropic model does 
not have a rigid 111 interface, although it does possess gapless excitations, 
as is well-known from spinwave theory. In the limit q —>■ 0, the Ising limit, 
the bound vanishes. This is to be expected, as the 111 interface contours of 
the Ising model are highly degenerate. 

In order to explain the role of the exponent 6{q, v) in ( |1.2[ ) we first need 
to discuss some properties of the interface states themselves. For < g < 1, 
the model has a two-parameter family of pure ground states with an interface 
in the 111 direction. One parameter is an angle, playing the same role as 
the angles (px in the Ansatz ( |1.4|) for the excitations. The second parameter, 
which is relevant for the present discussion, corresponds to the mean position 
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of the interface in the lattice. If we think of spin up at any site as describing 
an empty site, and spin down as a site occupied by a particle, the third 
component of the spin becomes equivalent to the number of particles. In 
Section 2, ( |2.8|) , we will introduce the chemical potential /i to control the 
expected number of particles, alias the third component of the total spin. In 
the limit q ^ 0, the filling factor u of the interface has a simple interpretation: 
z/ = means that interface separates a region entirely filled with particles 
from a region that is empty. A non-zero u means that there is a partially 
filled plane in between the filled and the empty region, with filling factor 
z/. It turns out that the exponent i/), can be considered as a function 
of /i alone. For each value of /i G M, we get an interface state, and S is 
the distance of /i to the integers, i.e., 6{fi) = min(|/i — [/^J |, |1 — /i + [/ij |), 
where [/iJ is the integer part of /x. In general, the relation between fi and u 
depends nontrivially on q. But for all g, < g < 1, one has 6{q, 1/2) = and 
0) = 1/2. For further details on the interdependence of the parameters 
q,S,fi, and z/, we refer to Section 6.1. 

We believe that 0{1/ B?) is the true behavior of the low-lying excitations. 
There are indications in the physics literature that this should indeed be the 
case 0. Our rigorous bounds are obtained using the variational principle: 
If 'i/'o is a ground state of and ip is any other state that is linearly 
independent of ipQ, then 

7 — ^1-^0 < TT-To : l/,/,J,M|2 • (1-'^) 



VMM, 

The first factor in the RHS is the energy of the perturbed state ifj. The second 
factor is necessary to correct for the non-orthogonality of ip and the ground 
state. In general, one would need to consider the orthogonal complement of 
ip to the entire ground state subspace of H^- In the present case however, we 
know that for each eigenvalue of the third component of the total spin, 
there is exactly one ground state. As we will only consider perturbations 
that commute with J*^^\ it is sufficient to take the orthogonal complement 

of 1p to 1pQ. 

Our ansatz for is of the following form 



^= n 
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The energy of such a state can be written as follows 

(^Ali^Al^^ ^ p^^jl_eos(</),. -</),)]. (1.5) 



fl = l 



where the P^^y are probabilities determined by the interface ground state. 
Px^y can be interpreted as the probability that the bond (x, y) belongs to "the 
interface contour", i.e., one of the sites is occupied by an up spin and one by 
a down spin. These probabilities decay exponentially fast as a function of the 
distance to the expected location of the interface. In particular, this shows 
that the interface is rigid and that the problem of calculating its excitation 
energies is quasi two-dimensional. In fact, the next step in our proof makes 



this explicit. We consider excitations of the form ( |1.4|) with 



0. = '50(^), R>1 

where 5 is a suitable scale factor, is a smooth function with compact 
support in M^, and x± is the component of x G Z^, orthogonal to the 111 
direction. It is shown that the energy of such excitations satisfies the 
bound 



7i? 



< 



C{q) ||V0||i. 



In principle, is a map from to the circle, and as such could have non- 
trivial topology. As we will only be considering small perturbations, this will 
be of no relevance here. It is, therefore, natural to take for an eigenfunc- 
tion belonging to the smallest eigenvalue of —A on a circular domain with 
Dirichlet boundary conditions, which minimizes of the Rayleigh quotient on 
the RHS, i.e., the Bessel function Jq. This is different from the so-called 



superinstanton Ansatz of Patrascioiu and Seller in [12], where they use the 
fundamental solution of the Laplace equation, instead of an eigenf unction. 

All our results are for ground states that are eigenstates of the third 
component of the total spin, which is a conserved quantity, and for ther- 
modynamic limits of such states. We will call this the canonical ensemble. 
Our derivation, however, relies on an equivalence of ensembles result for the 
interface ground states of the XXZ model. The state of the "small" volume 
A^, immersed in the much larger volume A, is well approximated by a grand 
canonical state with suitable chemical potential (see Chapter 2 for the precise 
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definitions), which does not have a fixed magnetization. As expected, this 
equivalence of ensembles holds only for observables that commute with the 
third component of the total spin which are analogous to the gauge invariant 
observables in particle systems. This equivalence of ensembles result is non- 
trivial. Although we only give the proof in dimensions 3, it is straightforward 
to generalize the proof to all dimensions > 3. Equivalence of ensembles (in 
the above sense) does not hold for the one-dimensional model. This can be 
derived from the results in In two dimensions, our method without mod- 
ifications, yields the equivalence of ensembles for volumes that grow as ^J~L 
in the 11 direction and as L in the direction of the interface. With additional 
work one can obtain equivalence of ensembles result for standard sequences 
of increasing volumes. 

As another application of equivalence of ensembles we prove the existence 
of the thermodynamic limit of sequences canonical ground states with a given 
density, i.e., magnetization per site. 

Concerning the gap above diagonal interface states in dimensions other 
than three we can make the following comments. First of all, diagonal inter- 
face states exist in all dimensions [IH. In one dimension there is a spectral 
gap above the ground states 0. In two dimensions an upper bound of order 
1/i? was proved in ^j. The method of this paper can be used to obtain a 
bound of order 1/i?^ also in two dimensions. In all dimensions greater than 
three our method can be applied without change to obtain equivalence of 
ensembles, the existence of the thermodynamic limit and an upper bound of 
order 1/i?^ for the excitation energies. 

The paper is organized as follows. Chapter 2 introduces the model and 
the geometrical setting. Chapter 3 deals with the equivalence of ensembles 
result which is a main ingredient of our proofs. The bound on the excitation 
energy is a product of two factors as in ( |1.3| ). A bound on the first factor, 
called the energy bound, is derived in Section 4. The second factor requires an 
estimate for the inner product of the ground state with the perturbed state, 
which is derived in Section 5. In Section 6 we prove a number of results for 
the grand canonical ensemble in one dimension that we use in the paper. 

2 Interface states of the XXZ model 

Our magnet occupies a volume A which is a subset of Z^. Let 61,62, 63 denote 
the standard basis vectors in Z^. (See Figure 0.) We let l{x) denote the signed 
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distance from the origin: l{x) — + x'^ + x^, where x — {x^, x^, x^) e 1?. 
Then 

5(A) = {(xo, xi) : |xo - xil = 1, /(xi) = /(xo) + 1} (2.1) 

describes the set of oriented bonds in I?. The infinite stick is, by defini- 
tion, the set of vertices of the form 

... - 62 - 63, -63, 0, 61, 61 + 62, 61 + 62 + 63, 61 + 62 + 63 + 61, . . . 

For any even integer L, the finite stick Eq of length L + 1 is then given by 

Eo = {x e E^ I -L/2 < l{x) < L/2} . 

We will take for A is a cylindrical region whose axis points in the 111 direction, 
where by cylindrical we mean that A can be obtained from a subset V of the 
l{x) =0 plane, which we will call the base, by adding to all vertices a; e F 
the finite stick Eq: 

A = {x + y I x e F, y e Eo} 

The equation l{x) — c, for any constant c, defines a cross-section of A, which 
contains exactly A = |F| vertices. Hence, |A| — (L + 1)A. We refer to these 
cross-sections as planes. 
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Figure 2: The projection onto the 111 plane of a cylindrical volume A with 
triangular base. The shading of the vertices depends on the value of l{x) modulo 3. 
The orientation of the bonds is indicated by arrows. Observe that each site has an 
equal number of incoming and outgoing bonds. 



As an example, the projection onto the plane l{x) = 0, of the vertices of 
A with triangular base is shown in Figure 0, with different shades depending 
on the value of l{x) modulo 3. The orientation of the bonds is indicated 
by arrows, and one may observe that each site on the interior of A has an 
equal number of incoming and outgoing bonds. By construction, A can be 
decomposed into one-dimensional sticks running parallel to the cylindrical 
axis, which we will generically call S. (See Figure |^.) One should observe 
that E is comprised entirely of nearest-neighbor pairs so that every site on E 
is connected to every other site by a sequence of bonds. This will allows us 
to exploit the well-known properties of the one- dimensional Heisenberg XXZ 
model to describe E. The Hamiltonian for the spin-| ferromagnetic XXZ 
Heisenberg model is given by 

Ha= E (2-2) 

{xo,xi)eB(A) 

where 

hi,., = -A~\sSsi':> + SS;S^^^)-SSsi'^ + U^-A{A^^^^^^ (2.3) 
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Figure 3: The bonds connecting the vertices of a stick S form a one-dimensional 
subsystem. 



and A > 1 is the "anisotropic coupling", A{A) 



'1 - 1/A2, and q, 



< g < 1, is the solution of A = |(g + g ^) The matrices S-, 
are the Pauli spin matrices acting on the site x, 



(a) 



1,2,3) 



1/2 
1/2 



-i/2 
i/2 



5(3) 



1/2 
-1/2 



(2.4) 

The terms containing A{A) cancel on all sites except at the top and bottom 
plane of the cylinder. The usefulness of the nearest-neighbor Hamiltonian 
stems from the fact that its action on any bond is given by 



h'\w--—zTiQ\w-\n)), 



q + q 



—,m-<i''\n)) 



q + q- 

In other words, h'^ is the orthogonal projection on the unit vector 

1 



-Mm)-\n)). 



(2.5) 
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There is a (|A| + l)-fold degeneracy in the ground states with a unique ground 
state for each value of total third component of the spin XlxeA ^^'^^ ■ "^^^ basis 
vectors of the Hilbert space (C^)®'^' can be labeled with particle configura- 
tions a = {a{x)}xeA, where a{x) is or 1, corresponding to |t) and ||), 
respectively. We write N for the operator defined by 

and let A{A,n) denote the collection of all configurations with N{a) = n. 
Following Ijl]] the ground states are given by 

MKn)= (8)?'^^^°^^M«(^)), (2.6) 

Note that the weights of a are invariant under any permutation of the sites 
for which planes are invariant. These states describe an interface located, on 
the average, in the plane determined by {L/2 + lx)A = 0. 
We denote ||'0o(A, n)||^ by Z{K,n). This quantity is given by 

Z(A,n)= -Q ^2/(x)a(x) (2.7) 

aG^(A,n) xeA 

We will treat ^(A, n) as a canonical partition function. It will be useful to 
consider, also, its grand canonical analogue: 

L 

z^^(A, ii) = Y ^(^' ^)?"'''" = + q^^^'^"^''"^)- (2.8) 

?i=0 xsA 

Then it is easily seen that Z'-^'^(A, yu) is the squared- norm of the grand canon- 
ical vector defined by 

|A| 

^«^(A,/x) = Y,q--^UKn) = (g)(|T) + g'(^)-^|i)). (2.9) 

n=0 xeA 

Due to the product structure, the thermodynamic limit is simply given by 
for all local observables X. 
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3 Equivalence of Ensembles 

A key step in our argument is the development of an equivalence of ensem- 
bles. Specifically, we will show that for a gauge-invariant local observable the 
canonical expectation is close to the grand canonical expectation for some 
suitably chosen chemical potential fi. Here n only depends on the total spin 
of the canonical ensemble, not on the form of the observable. From this, 
naturally follows a thermodynamic limit for gauge-invariant observables. We 
begin with activity bounds that show that the ratio of two canonical parti- 
tion functions with different particle numbers is approximately exponential 
in the difference of the particle numbers, i.e.. 



Z{A,n-k) ^ Z{A,n)q- 



-2kiJ, 



for I A; I -C n. More precisely, we have the following lemma. 

Lemma 3.1 (Activity bounds) For every volume k, |A| — {L ^ 1)A, the 

ratio of canonical partition functions for different number of particles can be 
bounded from above and below by activity bounds as follows. Let Aq be any 
constant. Suppose n, < n < A{L + 1), and fi are such that 



n-A{n)%l<\A,A'l\ (3.1) 



Then, for every k satisfying 



one has the bounds 



\k\ < IaoA^/^, (3.2) 



2 



and 



Z{A,n-k) - "^^^0'^^^ 

^) > C(Ao,A)-l/'['^-2<^)iS+2M«'^'-^l/('"^'), 

^(A, n — /c) 



(3.3) 



(3.4) 



where a = 2| lng|, 



L/2 
l=-L/2 



^,-^.;o^Osh\l{l-^l)y 
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and ^ 

C(Ao,A) = i±^. (3.5) 

Moreover, if fx is the solution of ^ — (N)^^^ = 0, then, also using the bounds 
for 0"^ given in l \6.1^ ), we obtain 

fc^(i-q^) Z(A ri) 2fc^(i-9^) 

C{Ao/2,Ar\-^<^ < q-^^^—^^—L- < C{A,/2,A)q-^^. (3.6) 

Z(A, n — k) 

Alternatively, if solves — (N)^*^ = 0, then we obtain 

fc^(i-g^) Z(A,n) 2fc2(i 

CiAo/2, A)"igMTT7^ < ^ ' \ < C(Ao/2, A)g^^. (3.7) 

Z(A, n — k) 

Proof: This can be obtained as follows. Let consider the grand canonical 
probability 

p(/.,n)=g-^-H^(^. (3.8) 

with 

Z{n) = J2 ^"^"^ (3-9) 

where Sj is the i-th one dimensional stick that we are decomposing our vol- 
ume in, and where Zccif^) is the grand-canonical partition function. Clearly, 
we have 

Zin) = • (3-10) 

n:|n|=n 



Define 



and we have 



P{f^,n) = p{f^,n) , (3.11) 

n:|n|=n 



(3.12) 

Z[n — k) p[fi,n — k) 

The idea now is to make use of the local central limit theorem for the proba- 
bility distribution of the occupation number in the i-th stick (see Theorem 
XVI.4.3.). Let = X^zGSi For any integer N, consider, the probability 

P^^i^i = ni, = un) =p(/i, n). (3.13) 
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Due to the factorization property of n), the ^'s are independent identi- 
cally distributed random variables. For centered i.i.d. random variables Xi 
with variance a^, the local central limit theorem guarantees that the random 
variable 

1 ^ 

^^ = ^TfJ2^- ■ (3.14) 
is close to a Gaussian in the sense that the quantity 

N 

Pn{x) := Prob(^X„ = x) (3.15) 

n=l 

fulfills the bounds 

1 :c2/c\ 1 rc2/c\ 

1 _ < p^(a;) < — 1 + —= (3.16) 



where c is the constant 

max(|x|, |x — fc|) 



(3.17) 



By applying ( p.l6| ) to the centered quantity Xn = in — {in), we obtain the 
following bounds on the ratio of probabilities: 

~ Pm{x -k) ~ 

where 

In terms of the non-centered variables ii we have 

p{fi,n) = P^(n-A(N)g3 (3.20) 

where (N)^^ is the average number of particles of a ID stick E, in the grand 
canonical ensemble with chemical potential fi. From this and the hypotheses 
(lO), (O), we obtain 
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Note that in case fj. is chosen so that (N)^*^ = n/A or (N)^*^ = {n — k)/A 
then we can replace c by c/2, with the result that C{Ao, A) may be replaced 
by 

C{Ao/2,A) = -±^, 

as well. 

Also, from and (|3J8| ), we have 

C(v4o,A)-ie ^^T^ < < CiAo,A)e 2.^a . 

p{lj,, n — k) 



(3.22) 

< C(Ao,A)e-'=[25-2<^>§.M+2-V-l]/2.^ ^ (3.23) 



Using ( |3.12| ) (and observing that q^^^ = e "'^), we have 

Z{n) 
Z{n - k) 

and 



Zin) 



GC I 9„„2,,_ fc 1 /9^2 



,^ > C(Ao,A)e-'=^'^-'<^^^.M+^"-^'^-tJ/^-\ (3.24) 
Z[n — k) 

Changing to base q then leads to equations ( p.3|) and (|3.4|) of the theorem. 
By the derivation of Section f).2\ we have the bounds on the variance for the 
number of particles in a ID stick: 

4 1 — I — 

In conjunction with the remark about replacing C{Aq, A) by C{Aq/2, A), this 
gives equations ( p.6[ ) and (|3.7|) . g 

As an application of this lemma, let us consider the case where n is 
replaced by p|A| — no, k is replaced by p|Ao| — no and A is replaced by 
Aq := A \ Aq. This means that in the lemma A is replaced by A — Aq, and 
(n - k)/A is replaced by p(|A| - |Ao|)/(A - Aq) = p{L + 1). Then, direct 
substitution shows 

Z(A^,p|A|-no) 

^(AS,P|ASI) 

< C{Ao/2, A - Ao) g-2fc.e-M2p(L+i)-2(N)gc+^]/2.2^ ^3^26) 
Z(A^,p|A|-no) 
^(AS,P|A§I) 
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where we have retained fc, for the moment. If, further, we choose /i so that 
{^)t.4i = p{L+l), which is always possible (see Section |0|) , then, by equation 
(|3.7|), we have 

^2.. ^(Ao,P|A| - ^o) ^ ^(^^/2, A - Ao) e-^^, (3.28) 
^2,. ^(Ao,P|A| -^o) ^ ^(^^/2,A_A„)-ie-^(:^^. (3.29) 

Using our bounds for cr^, we have 

^2,. ^(A^,P|A| - ^o) < ^^^^/2, A - Ao) e-^(^'^^ft-o), (3.30) 

''"'^yfrar S ^ - Ao)-' e-SfS . (3.31) 

By our choice of /x, conditions ( p.l|) and ( p.2|) are satisfied as long as the 
order of L does not exceed the order of {A — AqY^"^. This estimate will be of 
use in the next theorem. 

Let denote the operator- norm of X restricted to the subspace of 

ground states. For observables X, localized in A and commuting with J^^\ 
\\X\\gs is also given by 

= sup |(X)A,n|- 

0<ra<|A| 



Theorem 3.2 (Equivalence of Ensembles) Consider two cylindrical vol- 
umes A and Aq, Aq C A, of the type defined in Section 2 (in particular 
\A\ = A[L + 1), |Ao| = Ao{L + 1)), and fix a total number of particles n\. 
Define p = nA/|A|. Suppose X is a local observable in the volume Aq, which 
commutes with J^^^ := S-x^ . Then we have 

IWA,n- Wg^J < 4Xhs, (3.32) 

where 

_ \n\A - Ao) + 2(1 + a')Al + 4 AAp 



2{A-Ao) q\A-Aoy/^-2A^' 
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a = 2| lng|, and the chemical potential fi is determined by the equation 

(N)gS = p(L + l). (3.34) 



In particular, for p = 1/2 the calculations of Section \6. 1\ will show that 
/i = 0. 



Corollary 3.3 (Existence of the Thermodynamic limit) 

(i) Suppose we have a sequence of pairs {Ak,nk) with cylindrical volumes 
and Afe in such a way that the length does not grow faster than the linear 

size of the base. Let Hk solve (N)9'p = rik- Then the convergence jik ^ IJ' 
guarantees the convergence, of {■)A^,nk {■)^3^, for all local observables X 
commuting with J^^^ : 

^ {X)^f^ (3.35) 

(a) Moreover, for any choice of n, we may find a sequence of pairs {Ak,nk) 
such that 

(X)a„„, ^ {X)^f^. (3.36) 

Proof: (Proof of Corollary) It follows from the monotonicity of (N)^*^ 
proved in Section ^A\ , that the equation 

(N)?^:,, = Uk (3.37) 

always has a unique solution for fi^- Then, (i) follows immediately from 
the inequality ( p.32| ), once we observe that e \ as A in the sense 

prescribed in the corollary. 

For (ii), take A^, with base Ak, and Uk such that 



where [x\ denotes the largest integer < x. Then, pk solving ( |3.37|) , is easily 
seen to converge to p, and ( |3.36| ) follows from (i). _ 

The interpretation of the condition jj^k ^ p itl (i) of the Corollary is that, 
not only does nk/\Ak\ converge to p = 1/2, but, more precisely 

nk = p\Ak\ + uAk + oiAk) . 
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The term proportional to \Ak\ guarantees that the interface is in the center 
of the volume, the second term fixes its filling factor. 

Proof: (Proof of Theorem |3.2| ) Let fi be determined by ( p.34[ ), and define 
H as follows: 

„ ^ Z(A,nA)g-^^^l^°l 

^ Z(AS,p|Ag|)Z^^(Ao,/i) ^ • ^ 

where A^ := A\Ao. We will obtain the equivalence of ensembles by combining 
two facts. The first is that H is approximately equal to 1, and the second is 
an estimate showing that 



\{X)A,n.^-{X)f;j<e\\X 



\gs 



But first, let us recall the definitions of the expectation of an observable X: 

^ (^^^(A,p)|X|^^^(A,/.)) 

Since X is an observable localized in Aq, we note that {X)1'^^ = {X)'^^^^. 
Moreover, we may decompose the grand canonical state into a superposition 
of canonical states: 

|Ao| 

^^^(Ao, fi) = Yl ^"'^""^^(Ao, no). (3.41) 

reo=0 

Since X commutes with J^^\ it does not have off-diagonal matrix elements 
between these canonical states with all different values of the total spin. 
Therefore, 

|Ao| 

(X)£ = Z«^(A,/i)-i Yl 9"''^""^(Ao,no)(X)A„,.,. (3.42) 

no=0 

Note also, that since we have a decomposition 

|Ao| 

7A(A,n)= 5^^(A\A o,n - no) ip{Ao,no), (3.43) 

"0=0 
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and using the previously described properties, we have 



|Ao| 



^ Z{A \Ao,n- no)Z{Ao, hq) 



no=0 



Z(A,n) 

|Aol 



(3.44) 



no=0 

^^ Z(Ag,n-no)Zg^(Ao,/i) 
q-2f^noZ[A,n) 



(3.45) 



This differs from the definition of only by the final factor, which is a 

ratio of partition functions hence amenable to our activity bounds. 
In fact, we have 



GC 



|Ao| 

Z«^(Ao,/x)-^ J2 WA„,„„Z(Ao,no) 

"0=0 

2fi{no-(no)) ^(Ao,n - Uq) ^ 



X 



X 



(3.46) 



where {uq) = (N)^^^^, which equals p|Ao| for our choice of /i. Thus we obtain 



|(X)a,„S - (X)SI < \\X\U\9\)f^,,, where 



„ ^ 2^(no-(no» Z{A^o,n-no) _ 

^ ^ ^(a§,Lp|Ao|J) ■ 



(3.47) 



Now we use the activity bounds ( p.30| ) and ( p.31|) , but replacing k by its 
actual value, (no) — uq. We arrive at the bounds 



(i-g^)(("o>-"o)^ 



9 < 91 ■■= C(Ao/2, A - Ao)e 2(i+,^)(a-ao) _ 



where 



9 > 92-=C{Ao/2,A-Ao)-'e ^^^(^-^o) 

1 + 



C{Ao/2,A-Ao) 



2tT2(yl_ylo)l/2 



1 



(3.48) 
(3.49) 

(3.50) 



2ct2(A-Ao)1/2 

Therefore, < max(|5fi|, 15-21) < \9i \ + \92\- 

We now use the triangle inequality and the fact that the exponent is negative 



19 



to obtain: 



\9i\< 



(i-';^)«''o>-"o)^ 



I — e 2{l+5^)(yl-Ao) 



\l-C{Ao/2,A-Ac 



(3.51) 



so that 

(l-9^)((no)-no)^ 

(|^7i|)ao,m < (1 - e" ^(^-^^^)(^--o) + C(Ao/2, A - Ao) - 1. (3.52) 
Similarly, 

(|^72|)ao,, < (1 - e"''' W-V' + 1 - C(Ao/2, A - Ao)-^ (3.53) 



We will use the Chebyshev inequality to control the expectation term in 
( |3.52| ). Specifically, for any B > 0, 

(1 _ e 2(i+,^)(A-Ao) )GC^ ^ Prob(2|no - {uq) \>2B) + 1 - e 2(i+9^)(^-^o) 

In Section we show that (g~^'"°"^"''^')Ao^^ < 2(2g~2)Ao_ Qne choice for B 
is a~^[ln(74 — Aq) + ln(2g~^)]. This gives the bound 

(-o-_(^o»\ g^ 2 + ^^5^ [2(1 + a2)A2 + \n\A - A, 



(I - q ^-^0 < 

~ A — A, 



2 + {l + a')Al + lln\A-Ao) 

- A^o ^'-''^ 

=: Ci{A,Ao,q) 

2iA-Ao) 



The leading order term in the bound is '^/^"^^"^"^ for fixed q, strictly between 
and 1. Also, let 



which is greater than both C{Ao/2, A - Aq) - 1 and 1 - C{Ao/2, A - Ao)"^ 
Then |(/)A,nS-(/)£| < {C, + C2)\\X\\g,. In particular, |(I)a,„S- (]I)g| < 
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(Ci+C2)||3||gs, which is to say that |H — 1| < C1+C2. Then, using the triangle 
inequahty, we have 

|(X)a,.-(X)^,5| < |1-S|-|(X)a,„| + |(X)a,.H-(X)S| 
< 2(Ci + C2)||X||g,. 

So, defining e = 2Ci(g, A, Aq, n) + 2C2(g, A, Aq), the theorem is proved. _ 

Note that the restriction to observables X that commute with the third 
component of the total spin J^^^ is necessary. E.g., the expectation of 
obviously vanishes in any canonical state, while it is easy to see, by direct 
computation, that it does not vanish in the grand canonical states. This is 
entirely analogous to the restriction to gauge invariant observables in particle 
systems. 



4 Bound on the energy 

In this section we will estimate the energy of a class of perturbations of the 
ground state ipo given in (TB). Let A and be two cylindrical volumes 
as described in Section 2, C A. E.g., A/j and A, may have triangular 
cross-sections (see Figure |I]). We will generally assume that the radius R of 
An is much less than that of A. We consider ip of the form 

^(A,n,0)= J2 0e^'^(^)"(^)g'(^)"(^) , (4.1) 



where supp((/)) C A/j. 

We will also suppose that 



^m.m) (4.2) 



where is a smooth functions of its variables and iS is a parameter, which we 
will eventually take to zero independent of R. The coordinates y^,tildey'^, 
are defined by 

r> 1 2 S 2 

and are to be viewed as rescaled coordinates for x along the plane perpen- 
dicular to the 111 axis. 
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There are two points to our assumptions on 0: First, that (j) is independent 
of the 111 component of x. Second, that cj) is associated to a scale- invariant 
phase (f) by (j){x) = R~^(f){x/R). Ultimately, the constant S will vanish. The 
leading term in our estimate of the gap is independent of S as long as S -C 1. 

Let be the projection of Aji onto the plane l{x) = 0, Aji = \Tji\, flji 
be the convex hull of F/j, and = {a; G : Rx G ^r}, the rescaled region, 
and let m(f2) be the area of Cl (for the standard Lebesgue measure on M^). 

We will also use the following notation: dy(j) and are the first- and 
second-derivative tensors of 0, and by the norm of a tensor we mean the 
maximum of the L°° norms of the components. 

Then we have the following theorem. 

Theorem 4.1 (Bound on ^^^ff\J^'^ ) Considering a perturbed state as in 



the energy is bounded by 



where 



I < \J1 < 2l±^ (^^^^^^ + Snum 1 (4 4) 



£num= ^^llf \\dl<i)\\L--\\dy(i)\\L-- (4.5) 



i?5 

is a correction to the main term which becomes negligible as R —>■ oo 



Proof: We begin by calculating how a two-site hamiltonian h\ acts on the 



perturbed state. 

We consider the decomposition of our lattice into the relevant bond h = 
{xq, Xi) and everything else A \ 6. Thus 

/i^ = lA\b®IU(eb|, (4.6) 
where C,b is the unit vector from (p.5|) on the pair b, and 



^(A,n) = ^^(A\6, n — Ub) ^ ipibjUb). (4.7) 



nh=0 



Here ■ip{b, rib) is as would be defined by (f4.1|) , but with A replaced by b and n 
replaced by Ub- For example i){b, 1) = g'(^o)ei<^(xo) |||^ + ^^x^) ^i<t>{x^) |||^^ 
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is orthogonal to il){b, 0) and ■0(6, 1), since lies in the sector of total spin 
1. And 

Now it is straightforward to see 

{^{A,n)\hl\^{A,n)) (4.9) 

= mA\b,n-i)rmmb,i))\' 

= < 2 n2 ^(A,^)^^(&)(l-cos[0(a;i)-(/)(a;o)]), (4.10) 



where we have defined 



Then we may write 



^ff^=f^T^ ^ P^(6)(l-cos[0(xO-0(xo)]). (4.12) 
^ ' ^ ly "T- y ; 6es(A) 

Actually, P''(6) depends on b only through Z(xo). So from here on, we'll write 
it as P'^{l{xo)), and observe the following: 

/ / I ^(9) I / \ 9 ^/2-l 3 

(4.13) 

where = {x e Ar : l{x) = I}. 

Let us estimate the term ^^gpi ^^^^(1 — cos[0(x + ej) — 0(x)]). We have 
an inequality 

1 - cos[(P{x + ej) - (P{x)] < h(j){x + ej) - (j){x)]^ (4.14) 
(which is actually an equality in the limit — > oo for our ansatz). Also, 

Y],<l>{x + e,) - <\>{x)f ^ = l^lV^^r (4.15) 



1=1 
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In fact, using the inequahty 

m+v)-mf-[c-^ymf\<\\di 



\d.r, 



(4.16) 



one may conclude that the error in (|4.15| ) is bounded by ^||<9|0||loo \\dy(j)\\ 
Incorporating this estimate into the inequality of ( [4.14| ), we have 



— cos 



{x + e,) - <P{x)]) < 



2i?2 



(4.17) 



Finally, as i? — oo, the sum over each becomes increasingly well-approximated 
by the integral over we is proved in Lemma [4.2| immediately following 
this proof. The lemma gives us a bound 



52 IF 



R\ 



xeri 



i?2 



m(i7) Jh 



IV,- 



(4.18) 

where is the Laplacian and p = -^2/3 is the maximum radius for the 
Voronoi domain. (Note that by its definition, as the trace of the second- 
derivative tensor, the Laplacian enjoys the bounds 



(4.19) 



which may be combined with error term in ( ^4.17| ).) Combining ( [4.1^ ) and 
( 4.19| ) gives us the theorem, modulo the term X1z^-l/2 -^''(0) which we 



derive the necessary in Lemma 4.3 



Lemma 4.2 Suppose T is a region in a regular lattice. For each x G F, let 

be the Voronoi domain of x with respect to the whole lattice, and let 
he the union of all the individual domains Vt^ . If f is a smooth function on 
VLy, then 

<p||V,/|U^(np) (4.20) 
where p is the maximum radius of a Voronoi domain. 



|r| m(fir) V 
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Proof: For each x G F, 

^1 



1 



This clearly leads to the bound 



1 f f d 
Jn, Jo dt 

Vyf{x + t{y - a;)) ■ (y - x) dt dy. 



<p(n,.)l|V,/|U^(Q.). (4.21) 



From this, the lemma follows easily. b 
Now, we will derive the necessary bound on 

L/2-1 
l=~L/2 

We will rely on bounds for similar quantities in the one-dimensional model 
proved in 0. 

Lemma 4.3 (Bound on X1/^^l/2 ^^(0) 

L/2-1 2 

E ni)<2j^. (4.22) 

l=-L/2 ^ 

Proof: Recall 

Z[A, n) 

The ratio of partition functions in the equation above is clear: It is the 
probability of finding one particle shared by the sites of 6, and n — 1 particles 
shared by the sites of A \ 6, conditioned on finding n total particles on A. We 
consider the operator 

n = iiA\b®(in)b(nib + iiT)b(iTib)- 
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Then 

ZiA\b,n-l)Z{b,l) ,y, 

and 

L/2-1 / L/2-1 \ 

E ^'(0 = ( E ) • (4-25) 

l=-L/2 \l=-L/2 I 

where = (xo,a;i), where /(xq) = /, and {xq^Xi) is a bond in the stick 
containing the origin, which we denote by Eq. The restriction of the state in 
A with n spins down is of the form 

L+l 



k=0 



where X is any observable commuting with J*^^) = XlaieSo ' ^-S-; 
y6(;), and the Ck are non-negative numbers summing up to one. We will 

now derive an upper bound for (^;^^^/2 ^)so5 ^^^^ independent of the 
coefficients c^. We start from 

{Yih,,k < Prob,(5f ) =T, Sil\ =i) + Prob,(5f ) =i, =T) (4.26) 

where Prob^ denotes the probability in the ground state with k spins down 
for a one-dimensional system on [— L/2,L/2], the sites of which we label by 
/. Each term in the RHS of ( [4.26| ) can be estimate as follows. 

Probfc(Sf ) =T, =i) < min (Probfc(Sf ^ =T), Prohk{Sil\ =i)) (4.27) 

Theorem 7.1 of P] gives the following bounds 

Probfc(5g\ =i) < g2(/-(fe+l-L/2) if l>k+l-L/2 

Probfc(5f ^ =T) < g2(fe+i-i/2-0 if l<k+l-L/2 
Combining these inequalities and summing over / yields 

l=-L/2 ^ 

for all /c = 0, . . . , L -|- 1. Together with ( |4.25| ) this concludes the proof, g 
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5 Bound for the denominator 

Note that '?/'(A, n) = T{(j))ipo{A, n), where T(0) is the unitary operator defined 
by, 

n0) = (8)(lT)(Tl+e^*(^Mi)ai)- (5.1) 

In particular, \\T{(j))ipQ{A,n)\\'^ = ||'?/'(A, n)|p = Z{A,n). For convenience, 
we will sometimes omit the arguments A and n from the notation. In this 
section we will consider the half-filled system, i.e, p = n/\A\ = 1/2. This 
corresponds to /i = 0. 

Theorem 5.1 (Bound on |^^||^;^| ) Considering a perturbed state in the 
volume Aq defined by we have that canonical and grand- canonical ex- 

pectations of the perturbed state are arbitrarily close for large volumes A in 
the sense: 



(^ol^o) 



GC 



\n\A- Ao) + 2(1 + a^)Al + A AAq 



< — ^ ' . \\ " + 



2{A-Ao) g2^i/2_2Ao' 

(5.2) 

Moreover, with the ansatz defined by (§TJj, the grand canonical expectation 
is bounded as 

ln|(T(0))Sr< (5.3) 



< _o2'5(m) 



where S{p) is the distance of fi from its closest integer neighbor. (Recall that 
we have defined the L°°-norm of a tensor to be the L°°-norm of its maximum 
component.) 



Proof: The proof of equation ( p. 21 ) is a direct consequence of the equivalence 
of ensembles because, since T{(f)) is a unitary operator, ||T((/))|| = 1. Let us 
now consider the proof of equation ( |5.3| ). 

We wish to bound the denominator from below; i.e. to demonstrate 
that 1 — |(T(0))A,nP is not too small. This is tantamount to showing that 
|(T(0))A,nP is not too close to 1. Furthermore, we know this quantity lies 
between and 1. We estimate the actual canonical average with the grand 
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canonical average, and take the logarithm in order to exploit the factorization 
properties of the grand canonical ensemble. First, we note 



GC\ 



1 _|_ g«</'(^)g2{«(x)-/i) 
11 1 + nm^)'i^) 

xeAo 



(5.4) 



Recall the definition a = — 21ng. This allows us a more convenient form in 
place of (|]4D 



I _|_ g«0{3;)g2(/(a;)-/i) 



q 



2{l{x)-ii) 



n 

xeAo 



,2a{i(x)~i,) + 2cos0(x)e"('(^)-^') + 1 
Q2a{l{x)-tM) _|_ 2e"-(.K^)-t^) + 1 



(5.5) 



= JJ M - ^(1 - tanh2[a(/(x) - /i)/2])(l - cos0(x)) j . 
We partition the product over planes and estimate the logarithm, thus: 

ln|(T(0))£|' = ln(n 1-^(1 -ta^h'W^(^)-/^)/2])(l-cos0(x))) 

\xeAo / 

< ^ (l-tanh2[a(/(x) -/i)/2])(l-cos0(x)) 

xeA-o 

L/2 



l=-L/2 



We may approximate 1 — cos(0(x)) by |</)(a;)^, with an error no larger than 



^Il0lli°° which is the same as aiR? ^^is case 



In 



ZGc(Ao,/i, 0) 



ZGc(Ao,/i,0) 



^ L/2 

^ -9 E (l-tanh2[a(/-/x)/2]) 



l=-L/2 



E 



1,2 _ ^"^I^hI II I||4 



(5.6) 

We may approximate the sum over with an integral such that the error 
is bounded by — ^r^^H Vj;j0||Loo ||0||2,oo. We may bound the sum ^ 



L/2 

l=-L/2 



1- 
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tanh^[a(/ — /i)/2]) from below by its largest term (since all the terms are 
positive). The largest term occurs for that integer / which is closest to fi. 
Thus, defining 6{fi) = min(/i — \ — /i), we see 

J2 (l-tanh2[a(/-^)/2]) > l-tanh'[a6{fi)/2] = . ,(^) ,(^)., > q^'^'K 

l=-L/2 ^ ' 

(5.7) 

Using these bounds, we may continue the estimate of (|5.6|) . We arrive at 
ln|(T(0))^S|'< (5.8) 



< _g25(M) 



R\ 



^'(^) /^lly III Hill ^ 



Since |Vy0| < 2||9y0||;oo and since p = a/3/2, we have equation 



5.1 Bound on the Ratio 

We will now combine the results of the bound on the numerator and the 
bound on the denominator to get a true bound on the spectral gap. We 
first allow A in the appropriate fashion so that £ \ 0. Then we 

consider the case that S ^ 0, holding R fixed. This means that we consider 
a perturbation to the ground state which is very small. But since the ground 
state has energy zero, the energy of the perturbed state is entirely due to the 
small perturbation. In fact it is proportional to the size of the perturbation, 
and from this we obtain a linearized (with respect to amplitude of (p) bound: 
In fact we have, combining ( |1.3|) , (|4.4|) , and 



^ lg^2(l-5M) ||V50||2,(-^/m(fi) + |||9?0||oo||9^0||oo 

" (W)^ ■ ||0||J,(^/m(f2)-f ||a,0|U||0|U ^'^ 

Note that this bound is homogeneous with respect to the amplitude of (j), 
which is the result of our linearization. We observe that, whatever the form 
for 0, as long as it is smooth we have the same asymptotic behavior for the 
bound on the spectral gap. Namely 71 = 0(l/i?^). This said, it is certainly 
worthwhile to find a best bound, which we take up presently. 
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5.2 The Bessel Function Ansatz 

Let us write the leading-order term in the bound for the spectral gap: 

E(^) ^ SM. (5.10) 

In order to minimize the bound on the spectral gap, we will minimize the 
functional E{(j)) amongst all functions which possess two continuous deriva- 
tives and which vanish on the boundary of the rescaled perturbed region fl. 
(In order that the "small" phase match the external phase of 0, ±27r, . . . 
on dQ, it must be zero there. Thus = on dCl.) Therefore, we consider 
the first variation 

lim-|g(^ + r^')-gW]= " - 'jlJ// - (5.11) 

Setting the first variation to zero for all test functions 0' leads to the eigen- 
value problem for Laplace's equation 

-VV = A0 inn, 

i) = o on an, ^ ■ ^ 

where A = E{(j)). 

We choose, for our domain, the unit disk. We seek the solution to equation 
( |5.12| ) which minimizes A, but with the restriction that must possess two 



continuous derivatives. So the fundamental solution, which is the logarithm, 
is disallowed (and, in fact, has higher energy). We seek the first eigenstate 
of the Laplacian above the ground state. This is a classic problem, found in 
any elementary PDE text, with the Bessel Function for the solution: 

where r = \y\, Jq is the zeroth Bessel function, and zq ~ 2.406 is its first 
zero. Now, using this choice for and the bounds (|5.9|), we obtain 



IQq^a-SM) L56 + |(2.90)(L40) 
- (1 - g2)i?2 ■ 27-f (L40)(l) ' ^^'^^^ 



Thus, 

lOOg^d-SCM)) 



71 < (i_g2)^2 ^>70. (5.14) 
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6 Results from the ID grand canonical en- 
semble 

6.1 The mean number of particles in a stick 

Recall that E is a ID stick running parallel to the 111 axis. So, it is actually 
a ID spin chain. We wish to estimate the mean number of particles in E, for 
the grand canonical ensemble. This is 

L+1 



(N)gJ := Z«^(E, /.)-^ Yl nq-''''Z{E, n) (6.1) 

n=l 
L+1 

= Z«^(E, ix)-^ ^ ne"'^"Z(E, n). 

n=l 

where E is the interval { — f , — f + 1, . . . , f }. (Recall a — — 21ogg.) By a 
standard calculation, we have 

(N)S = ^|^log^^^(S'/^)- (6-2) 

On the other hand, the grand canonical partition function factorizes, as we 
have seen, so that 



(l^)gS = E = E i [l - to-h (f (' - /-))] ■ (6.3) 

l=-L/2 l=-L/2 



An examination of the graph of the function a; i— > 1 — tanh(x) reveals an 
approximate heaviside function, with support on the negative axis. We define 
the function 

{1 X < 0, 
1/2 x = 0, (6.4) 
X > 0. 

Then, as long as — L/2 < ji < L/2, we remark 

(N>is i }+ _E (i - \ (if - - - ")) . 

(6.5) 
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We make the definition 



L + 1 



(6.6) 



For /i in the range above one may determine (by combining the two tails in 
the series and estimating upwards by an integral) that 



a V 1 + exp( 



1 , /l + exp(-f(f -/x)) 
+ /i)) 



L 
2 V 2 
afL_ 
"2V2 



(6.7) 



Notice that in case /x = 0, there is no error at all in estimating Fl by Foo, 
and, furthermore, -Foo(O) = 0. It is clear that -Foo(/i) is periodic in /i with 
period 1, because it is a sum over the entire integer lattice, so it will suffice 
for us to consider /z in the range ]0, 1[. A straightforward calculation then 
yields 



Foo(/i) 



i i ^ 

-u H h > 

—fl H — tanh(aju) + 



1=1 

oo 



1 



I _|_ ^a{l-i^i) ]_ _|_ ^a{l+ti) 

sinh(a/i) 



Defining {^} 

FM -- 



^ — we have 
1 



cosh(a/i) + cosh(a/) 
sinh(a{/i}) 



{/x} + - tanh(a{/i}) + V bumw/xj; 
^'^^ 2 ^ cosh(a{^}) + cosh(a/) ^ ^ 



for all values of fi. 

Lemma 6.1 The function F^o defined in ( |6'.<^ has the following properties: 

i) Foo is periodic with period 1, i.e, Fooifi + 1) = Foo(/u), for a// G M. 

ii) Foo is odd about ^ = 1/2, i.e., F^oil — fi) = —Fooifi), for all /i G M. 
Hi) -1 < Foo(Ai) < 1, for all G M. 

iv) Foo(/i) = for /i G Z and G | + Z. I.e. the estimate (N)^ 
is exact for half-integer and integer filling. 



\GC 



Proof: The periodicity of F^o follows directly from its definition. To prove 
(ii), define F{fi) for < /i < 1 as 



Fif^) = E 



k=l 



1 



1 



X _|_ ga(«-At) X + e"('+^) 



1 



1 + e"'^ 



(6.9) 
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Then, 



oo 



1 1 



1 



X _|_ ga(«-l+^t) ]_ _|_ ga{i+l-^t) 



1 _(_ ea(l-/i) 



«=1 




X -|- ga{«+M) 1 + e"('~^) 



11 1 




And clearly the remainder term 




if /i G Z 



satisfies property (ii). For the bounds, we first restrict ourselves to fi E [0, 1]. 
For /i > 0, we note that (|6.8|) implies 



Then we use property ii) in combination with this bound to also get the 
upper bound for /i G [0, 1]. 



Due to the peridicity property i), the upper and lower bound are automati- 
cally extended to all real /i. The special values stated in iv) are straightfor- 
ward from ( |6.8|) and ( |6.9|) . g 

We can define the quantity S{fi) = min(|/i — [//J |, |1 — /i + [/^J |), where 
[/ij is the integer part of /i. In general, the relation between /i and u depends 
nontrivially on q and the function 6 can be thought as S{q, v). But for all q, 
< g < 1, one has 8{q, 1/2) = and 0) = 1/2. See Figure |. 

6.2 The variance of the number of particles in a stick 

In the same way as was done above for the mean, we can compute the variance 
of the number of particles in a stick in the grand canonical ensemble by using 
the standard formula 



a\^^^L) = (N^)gS - ((N)^^)^ = log ^), (6.10) 



FM > -M > -1. 



i"oo(l-/i) < 1 
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Figure 4: A plot of the functions F^{fi) and cr {^i), with q = e 
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which gives 



Define 



L/2 



l=-L/2 



cosh'(|(i-//)) 



Then, the speed of convergence of this hmit is bounded as follows: 



2q2{L/2-n) 



n=0 



(6.11) 



(6.12) 



(6.13) 



It is clear that cr^(/x) is a periodic function of /i with period 1. It is not hard 
to see that cr^(//, L) is C°° and attains its maximum in all integers and its 
minimum in the integers +1/2. It is easy to derive upper and lower bounds 
for (7^(/i, L). An upper bound is given by 

(j\ii,L)< e-K'-'^)l< e-«l'l<l + ^^ (6.14) 

l=-L/2 l=-L/2 ^ 
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Figure 5: A plot of the function d^i'jq) for four different values of q. 



and a lower bound can be obtained using the crude bound 2 cosh a; < 2e'^': 
.U^)>;E^-'°"'s f 1,V„ (6.15) 

n=l 

From ( |6.14|) and ( |6.15|) we see that the limit satisfies the bounds 

71 2 < ^ < 1 6.16 

for all real fi and where we have again used the relation e""^ = g^. 
For the afficionados, one can also show that 

lim(T^(u) = ji '.f^^rj, (6.17) 

The interpretation is simple. When yU G Z, the interface (kink) in the one- 
dimensional system is located at a lattice site, which is occupied by a particle 
with probability 1/2. Clearly, the variance of the particle number is them 
1/4. However, for yU ^ Z, the kink is centered at a position not belonging to 
the lattice and the state converges, as g J. 0, to a deterministic configuration 
with zero variance for the particle number. 
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6.3 Estimating (g2|N-(N)l)GC 
We begin with the obvious fact 

from which it follows that 

Now, we observe 

Since 

Z^^(S,/i)= JJ L/2(l + g2{'"M)) (6.21) 

equation ( |6.20| ) leads us to conclude 

^ -L „2(l/2+i-m) 

{rh.,= \Xl_,,„,,^, <2,^W^^>-'. (6.22) 

Similarly, 

1 _L ^-2{L/2+1+m) 

= £ 2,-»/^«+'->. (6.23) 



Using the results of section |6.1| , we then have 

^^2|N-(N)|^GC < 4^-l-|F,(;.)| < 4^-2_ (q_24) 

If we wish to calculate (q'^''^~^'^^')a^, where A is comprised of A sticks, then 
nothing changes except that each estimate is raised to the power A. Thus, 
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